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Abstract 

We consider the dynamics towards the initial singularity of Bianchi type IX vacuum 
and orthogonal perfect fluid models with a linear equation of state. Surprisingly few facts 
are known about the 'Mixmaster' dynamics of these models, while at the same time most of 
the commonly held beliefs are rather vague. In this paper, we use Mixmaster facts as a base 
to build an infrastructure that makes it possible to sharpen the main Mixmaster beliefs. 
We formulate explicit conjectures concerning (i) the past asymptotic states of type IX 
solutions and (ii) the relevance of the Mixmaster/Kasner map for generic past asymptotic 
dynamics. The evidence for the conjectures is based on a study of the stochastic properties 
of this map in conjunction with dynamical systems techniques. We use a dynamical systems 
formulation, since this approach has so far been the only successful path to obtain theorems, 
but we also make comparisons with the 'metric' and Hamiltonian 'billiard' approaches. 



'Electronic address: mark.heinzleOunivie . ac . at 
^Electronic address: claes.uggla@kau.se 



1 



1 INTRODUCTION 



2 



1 Introduction 

Today, Bianchi type IX enjoys an almost mythical status in general relativity and cosmology, 
which is due to two commonly held beliefs: (i) Type IX dynamics is believed to be essentially 
understood; (ii) Bianchi type IX is believed to be a role model that captures the generic features 
of generic spacelike singularities. However, we will illustrate in this paper that there are reasons 
to question these beliefs. 

The idea that type IX is essentially understood is a misconception. In actuality, surprisingly 
little is known, i.e., proved, about type IX asymptotic dynamics; at the same time there exist 
widely held, but rather vague, beliefs about Mixmaster dynamics, oscillations, and chaos, which 
are frequently mistaken to be facts. There is thus a need for clarification: What are the known 
facts and what is merely believed about type IX asymptotics? We will address this issue in two 
ways: On the one hand, we will discuss the main rigorous results on Mixmaster dynamics, the 
'Bianchi type IX attractor theorem', and its consequences; in particular, we will point out the 
limitations of these results. On the other hand, we will provide the infrastructure that makes it 
possible to sharpen commonly held beliefs; based on this framework we will formulate explicit 
refutable conjectures. 

Historically, Bianchi type IX vacuum and orthogonal perfect fluid models entered the scene 
in the late sixties through the work of Bclinskii, Khalatnikov and Lifshit^] [U [2] and Misner 
and Chitre [HI HI El E] - BKL attempted to understand the detailed nature of singularities and 
were led to the type IX models via a rather convoluted route, while Misner was interested in 
mechanisms that could explain why the Universe today is almost isotropic. BKL and Misner 
independently, by means of quite different methods, reached the conclusion that the temporal 
behavior of the type IX models towards the initial singularity can be described by sequences of 
anisotropic Kasner states, i.e., Bianchi type I vacuum solutions. These sequences are determined 
by a discrete map that leads to an oscillatory anisotropic behavior, which motivated Misner to 
refer to the type IX models as Mixmaster models [3J [4j . This discrete map, the Kasner map, was 
later shown to be associated with stochasticity and chaos (20 [9], a property that has generated 
considerable interest— and confusion, see, e.g., [IHl HH H3 UHl E31 Ull 113 UHl H3 HD] and 
references therein. A sobering thought: All claims about chaos in Einstein's equations rest 
on the (plausible) belief that the Kasner map actually describes the asymptotic dynamics of 
Einstein's equations; as will be discussed below, this is far from evident (despite being plausible) 
and has not been proved so far. 

More than a decade after BKL's and Misner's investigations a new development took place: 
Einstein's field equations in the spatially homogeneous (SH) case were reformulated in a manner 
that allowed one to apply powerful dynamical systems techniques [2TJ[22l|23]; gradually a picture 
of a hierarchy of invariant subsets emerged where monotone functions restricted the asymptotic 
dynamics to boundaries of boundaries, see [lOj and references therein. Based on work reviewed 
and developed in [TD] and by Rendall [53], Ringstrom eventually produced the first major proofs 
about asymptotic type IX dynamics [25l [26] . This achievement is remarkable, but it does not 
follow that all questions are settled. On the contrary, so far nothing is rigorously known, e.g., 
about dynamical chaotic properties (although there are good grounds for beliefs), nor has the 
role of type IX models in the context of generic singularities been established [27] [28l [29] [30] . 

The outline of the paper is as follows. In Section [2] we briefly describe the Hubble-normalized 
dynamical systems approach and establish the connection with the metric approach. For sim- 
plicity we restrict ourselves to the vacuum case and the so-called orthogonal perfect fluid case, 
i.e., the fluid flow is orthogonal w.r.t. the SH symmetry surfaces; furthermore, we assume a 
linear equation of state. In Section [3] we discuss the levels of the Bianchi type IX so-called 



We will refer to the authors and their work as BKL. 
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Lie contraction hierarchy of subsets, where we focus on the Bianchi type I and type II subsets. 
In Section [4] we present the results of the local analysis of the fixed points of the dynamical 
system and discuss the stable and unstable subsets of these points which are associated with 
non-generic asymptotically self-similar behavior. Section[5]is devoted to a study of the network 
of sequences of heteroclinic orbits (heteroclinic chains) that is induced by the dynamical system 
on the closure of the Bianchi type II vacuum boundary of the type IX state space (which we 
refer to as the Mixmaster attractor subset). These sequences of orbits are associated with the 
Mixmaster map, which in turn induces the Kasner map and thus the Kasner sequences. We 
analyze the properties of non-generic Kasner sequences and discuss the stochastic properties of 
generic sequences. In Section [6] we discuss the main 'Mixmaster facts': Ringstrom's 'Bianchi 
type IX attractor theorem' [25|, 126] , Theorem 16. 1[ and a number of consequences that follow 
from Theorem 16.11 and from the results on the Mixmaster/Kasner map. In addition, we in- 
troduce and discuss the concept of 'finite Mixmaster shadowing'. In the subsection 'Attractor 
beliefs' of Section [7| we formulate two conjectures that reflect commonly held beliefs about 
type IX asymptotic dynamics and list some open issues that are directly connected with these 
conjectures. In the subsection 'Stochastic beliefs' we address the open question of which role 
the Mixmaster/Kasner map and its stochastic properties actually play in type IX asymptotic 
dynamics. This culminates in the formulation, and discussion, of two 'stochastic' conjectures. 
In Section [5] we present the Hamiltonian billiard formulation, see 0(3] or [3T]; we demonstrate 
that this approach yields a 'dual' formulation of the asymptotic dynamics. We point out that 
the billiard approach is a formidable heuristic picture, but fails to turn beliefs into facts. We 
conclude in Section [9] with a discussion of the main themes of this paper. Throughout this paper 
we use units so that c = 1 and 8irG = 1, where c is the speed of light and G the gravitational 
constant. 



2 Basic equations 

We consider vacuum or orthogonal perfect fluid SH Bianchi type IX models (i.e., the fluid 4- 
velocity is assumed to be orthogonal to the SH symmetry surfaces) with a linear equation of 
state; we require the energy conditions (weak/strong/dominant) to hold, i.e., p > and 

-±<w<l, (1) 

where w = pj p, and where p and p are the energy density and pressure of the fluid, respectively. 
By ([T]) we exclude the special cases w — — 3 and w = 1, where the energy conditions are only 
marginally satisfied!! 

As is well known, see, e.g., [lOl [26] and references therein, for these models there exists a 
symmetry- adapted (co-)frame {w 1 , u> 2 , u> 3 }, 

du) 1 = —hi u> 2 A u> 3 , clu) 2 = -h 2 lj 3 A lj 1 , du 3 = -n 3 uj 1 A Cj 2 , (2a) 
with hi = 1, 722 = lj ^3 = 1; such that the type IX metric takes the form 

4 g = -dt® dt + g xx {t) w 1 + .922W & 2 <g> lj 2 + g 33 (t) uj 3 ® u> 3 . (2b) 

2 Note that the well-posedness of the Einstein equations (for solutions without symmetry) has been questioned 
in the case — l/3<iu<0, see I32| . The case w = 1 is known as the stiff fluid case, for which the speed of sound 
is equal to the speed of light. The asymptotic dynamics of stiff fluid solutions is simpler than the oscillatory 
behavior characterizing the models with range — i < w < 1, and well understood 26 33 . (In the terminology 
introduced below, the stiff fluid models are asymptotically self-similar.) We will therefore refrain from discussing 
the stiff fluid case in this paper. 
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Table 1 : The class A Bianchi types are characterized by different signs of the structure constants 
(n a , hp, n 7 ), where (a(3j) is any permutation of (123). In addition to the above representations 
there exist equivalent representations associated with an overall change of sign of the structure 
constants; e.g., another type IX representation is ( ). 



Hence, the type IX models naturally belong to the so-called class A Bianchi models, see Tablc[T] 
Let 

m{t) :=m —7== , n 2 {t) := n 2 j== , n 3 (t) := ti 3 — — , (3) 
Vdet g ydet g y/aet g 

where det g = 511522533- Furthermore, define 

= -trfe and a a = —k a p + |tr k £°g = diag(oi, cr 2 , <7 3 ) X/ ° a = °) ' ^ 

where k a p denotes the second fundamental form associated with ^ of the SH hypersurfaces 
t = const. The quantities 9 and a a p can be interpreted as the expansion and the shear, 
respectively, of the normal congruence of the SH hypersurfaces. In a cosmological context it 
is customary to replace 8 by the Hubble variable H = 9/3 = — trfc/3; this variable is related 
to changes of the spatial volume density according to d^/detg/dt = 3H^/det 5. Evidently, in 
Bianchi type IX (and type VIII) there is a one-to-one correspondence between the 'orthonormal 
frame variables' (H,a a ,n a ) (with ^ a a a = 0) and (gaPikap); in particular, the metric g a p is 
obtained from (ni, 712,713) via ©. (For the lower Bianchi types I- VIIo, some of the variables 
(711,712,713) are zero, cf. ([3]); in this case, the other frame variables, i.e., (H,a a ), are needed as 
well to reconstruct the metric; see [34] for a group theoretical approach.) 

In the Hubble-normalized dynamical systems approach we define dimensionless orthonormal 
frame variables according to 

(E a ,N a ) = (a a ,n a )/H, n = p/(3H 2 ). (5) 

In addition we introduce a new dimensionless time variable r according to dr/dt = H. Like the 
cosmological time t, the time r is directed towards the future; however, to make contact with 
the well established convention that uses a past-directed 'time' for the discrete Mixmaster map, 
see Section [SJ it will occasionally become necessary to use an inverse time r_ = — r instead of 
r itself. 

For all class A models except type IX the Gauss constraint guarantees that H remains positive 
if it is positive initially. In Bianchi type IX, however, it is known from a theorem by Lin and 
Wald [35 that all type IX vacuum and orthogonal perfect fluid models with w > first expand 
(H > 0), reach a point of maximum expansion (H = 0), and then recollapse (H < 0)H The 
variable transformation ([SJ breaks down at the point of maximum expansion in the type IX 
case; however, the variables (E a ,N a ) correctly describe the dynamics in the expanding phase, 
which we will focus on henceforth. 

3 In the locally rotationally symmetric case it has been proved that the range of w can be extended to to > — i, 
see |36| . There are good reasons to believe that the assumption of local rotational symmetry is superfluous, but 
this has not been established yet. 
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When the Einstein field equations are reformulated in terms of H and the Hubble-normalized 
variables (£ a ,AT a ) it follows for dimensional reasons that the equation for the single variable 
with dimension, H, 

H' = -(l + q)H, (6) 

decouples from the remaining dimensionless equations [ 10] ; here and henceforth a prime denotes 
the derivative d/dr. The equations for (T, a ,N a ) form the following coupled system [TO] : 



£' a = -(2-g)£ a - 3 S Q 
N' a = (q + 2£ Q ) N a 



where 



q = 2£ 2 + |(l + 3«))fi , 



(no sum over a) 



* 2 = m 



si), 



and 3 S a = ± [N a (2N a - N/3 - N 7 ) - (Np - N 7 ) 2 ] , (a/3 7 ) € {(123), (231), (312)} 



(7a) 
(7b) 

(8a) 
(8b) 



Note that the 'deceleration parameter' q is non- negative because of the assumption w > —1/3. 
Apart from the trivial constraint Si + £2 + £3 = 0, there exists the Gauss constraint 



12 



+ Ni + iV 3 2 - 2 (NiN 2 + N 2 N 3 + N 3 Nx) 



Q = 1 



(9) 



which is used to globally solve for f2 when ^ 0. Accordingly, the reduced state space is given 
as the space of all (£1, £2, £3) and (N\, N 2 , N3) such that £1 + £2 + £3 = and 



Nl + Nl + Nl - 2 {N X N 2 + N 2 N 5 + N 3 N X ) 



< 1 



®) 



which follows from © under the assumption that > (p > 0). It follows that the di- 
mensionless state space of the Bianchi type IX orthogonal perfect fluid models with a linear 
equation of state is 5-dimensional0 while the state space of the vacuum models (i.e., Q = 0) is 
4-dimensional. The same is true for Bianchi type VIII, while the state spaces of the remaining 
class A Bianchi models have less degrees of freedom; see Table [2] Once the dynamics in the 
dimensionless state space is understood, H is obtained from a quadrature by integrating ([5]), 
which allows one to reconstruct the metric. 

For all class A models except type IX the constraint ([9lJ implies that £ 2 < 1; in Bianchi 
type IX, however, £ 2 > 1 is possible. For type IX, define 



KiViiVaiVa) 2 / 3 



Employing (f5aj) and (O and using that 



Nl + iV 2 2 + iVf 



2A T 



A > 0. 



where equality holds iff Ni = N 2 = N3, we find 



£ 2 < 1 + A 



fi < 1 + A 



2-q> §(l-w)ft-2A 



(10) 

(11) 
(12) 



The function A is strictly monotonically increasing along orbits of Bianchi type IX. To see this 
we use ([7]) and compute 



A' = 2qA . 



A' 



= 



A'' 



9=0 



q=0 



= § [ 3 S 2 + 3 S 2 2 + 3 S 2 3 ] A , 



(13) 



4 It is common to globally solve Ej + £2 + S3 = by introducing new variables according to Si = — 22+ , 
22 = 2+ — \/32_, 23 = 2+ + \/32_, which yields 2 2 = 2^ + However, since this breaks the permutation 
symmetry of the three spatial axes (exhibited by type IX models), we choose to retain the variables 2i, 22, 23. 
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Table 2: The dimensionless state spaces associated with class A Bianchi models; here, (a/3j) 
is any permutation of (123). In addition to the above representations there exist equivalent 
representations associated with an overall change of sign of the variables (N±, N2, N3). The 
quantity D denotes the dimension of the state space (in the fluid case) ; the dimensionality of 
the state space in the vacuum cases is given by D — 1. 



where we note that 3 ,5 2 + 3 S 2 2 + 3 5 f 3 2 > because of the constraints. In combination with (Tj"2"|) 
it follows that S 2 and fl are bounded towards the past. 

The right hand side of the reduced system ([7]) consists of polynomials of the state space variables 
and is thus a regular dynamical system. Solutions of ([7]) of Bianchi types I- VIII are global in 
r, since ([S]) implies the bounds S 2 < 1 and fi < 1 which control the evolution of N a in (I7b[) . 
Solutions of ([7J of Bianchi type IX are global towards the past, since q + 2E Q is bounded from 
below; this follows from pa), which yields q > 2E 2 , so that q + 2E Q > -2 + ±(£ a + 2) 2 + 
g(£/3 — S 7 ) 2 . The decoupled equation for H, cf. ([B]), yields that -ff-^ooasr^— 00, because 
q is non-negative. Since q is bounded as t — > —00, the asymptotics of H can be bounded by 
exponential functions from above and below. It follows that the equation dt/dr = H^ 1 can be 
integrated to yield t as a function of r such that t — > as t - >— 00. 

In addition to (J7J) it is useful to also consider an auxiliary equation for the matter quantity f2, 

9! = [2q - (1 + 3w)]Q . (14) 

Making use of (O and ([5]) we conclude that for all orthogonal perfect fluid models with a linear 
equation of state we have p oc exp (— 3[1 -I- w]t), and hence p — > 00 as r — > —00, which yields 
a past singularity. The divergence of p can also be directly read off from the matter equation 
V a T ab = 0. 

3 The Bianchi type IX Lie contraction hierarchy 

The Bianchi type IX state space Bix is characterized by the conditions N± > 0, N 2 > 0, N 3 > 0. 
We write £?ix = Bn 1 n 2 n 3 - The notation is such that the subscript denotes the non-zero variables 
among {N\, iVjj, N3}. Setting one or more of these variables to zero (which corresponds to Lie 
contractions [37]) yields invariant boundary subsets which describe more special Bianchi types. 
Since the type IX models exhibit discrete symmetries associated with axes permutations, the 
contractions generate all possible representations of the more special (Lie contracted) Bianchi 
types (which are associated with such permutations): The Bianchi type VIIo subspace -Bviio 
is given by the disjoint union of three equivalent sets, -Bvn — Bn ± n 2 U Bn 2 n z U Bn 3 Nx-> where, 
e.g., Bn 1 n 2 denotes the type VIIo subset with N\ > 0, N 2 > and N 3 = 0; the Bianchi type II 
subspace B\\ by the union Bu = Bn ± UBn 2 U2?at 3 ; the Bianchi type I subspace B\ by B\ = Bj>. 
Note that the Bianchi type VIo subspace does not appear as a boundary subset of B\x- A 
Bianchi subset contraction diagram for type IX is given in Figure [TJ 
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Figure 1: Subset contraction diagram for Bianchi type IX. D denotes the dimension of the 
dimensionless state space for the various models with an orthogonal perfect fluid with linear 
equation of state; the associated vacuum subsets have one dimension less; see also Tabled The 
notation is such that the subscript of B+ denotes the non-zero variables, e.g. Bm 1 denotes the 
type II subset with N x > 0, N 2 = and N 3 = 0. 



Each set of the Lie contraction hierarchy is the union of an invariant vacuum subset, i.e., 
O = 0, and an invariant fluid subset, i.e., f2 > 0. To refer to a vacuum [fluid] subset of a 
Bianchi set B+ we use the notation Bl ac - [#*']. In this spirit, e.g., the type II subset decomposes 
as B Nl = B V N ^- u B R Ni . 

In the following we analyze the boundary subsets to the extent needed in order to understand 
the asymptotic type IX dynamics. 



The Bianchi type I subset 

The Bianchi type I subset is given by Ni = 0, N 2 = 0, N 3 = and ft = 1 - £ 2 > 0; since 
Ni, N%, N3 vanish, we denote this subset by B$, cf. Figure [TJ The vacuum subset consists of a 
circle of fixed points — the Kasner circle K°, which is characterized by £ 2 = 1. It is common 
to represent different points on K° in terms of the Kasner exponents p a , 

(Si, £2, S 3 ) = (3pi - l,3p 2 - i,3p 3 - 1) ; Pi +j?2 +i?s = 1 , pI+pI+pI = ±; (15) 

each fixed point on K° represents a Kasner solution (Kasner metric) with the corresponding 
exponents. The Kasner circle is divided into six equivalent sectors, denoted by permutations 
of the triple (123), where sector (a/37) is characterized by p a < pp < p 7 , see Figured] The 
boundaries of the sectors are six special points that are associated with solutions that are locally 
rotationally symmetric (LRS): Q a are given by (£ a ,£/3,£ 7 ) = (—2,1,1) or (p Q ,P/3,p 7 ) — 
(— i, |, |) and yield the three equivalent LRS solutions whose intrinsic geometry is non-flat; 
the Taub points T Q are given by (£ a ,£g,£ 7 ) = (2,-1,-1) or (p a >P/3>P 7 ) = (1,0,0) and 
correspond to the flat LRS solutions — the Taub representation of Minkowski spacetime. 

The quantity P1P2P3 (or, equivalently, £i£2£3 = 2 + TJp\p2Pz) is invariant under changes of 
the axes and thus naturally captures the 'physical essence' of a solution independent of the 
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Figure 2: The Kasner circle K° of fixed points divided into its six equivalent sectors and the 
LRS fixed points T Q and Q a . Sector (a/37) is defined by E Q < E^ < E 7 . The sectors are 
related to each other by permutations of the spatial axes. 



chosen frame; however, it is typically replaced by the Kasner parameter u through 

-it 2 (l + u) 2 

P1P2P3 = 77-. ; — 2vT ' where u € [l,oo] . (16) 

The Kasner parameter u parameterizes the Kasner exponents uniquely (up to the permutation 
symmetry); we have 

—u l + u u(l + u) 

1 + u + u z 1 + u + u A 1 + u + u z 

for sector (a/37) °f K°, where 116 (1, 00). Therefore, each point on sector (a/37) is represented 
by a unique value of u € (1, 00). At the boundary points of sector (a/37), which are Q Q and T 7 , 
the Kasner parameter is u = 1 and u = 00, respectively. Permuting (a/37) yields a physically 
equivalent state on a different sector; accordingly, each u S (l,oo) represents an equivalence 
class of six points on K°. In contrast, u = 1 describes the three points {Qi, Q2, Q3}; U — 00 
yields {Ti,T 2 ,T 3 }. 

While the Bianchi type I vacuum subset coincides with the Kasner circle K°, which is given 
by E 2 = 1, the Bianchi type I perfect fluid subset is the set 1 — O = E 2 < 1. From (fT4l) it is 
straightforward to deduce that there exists a central fixed point, the Friedmann fixed point F, 
given by E Q = Va, which corresponds to the isotropic Friedmann-Robertson- Walker (FRW) 
solution. Solutions with < E 2 < 1 are given by radial straight lines originating from K° and 
ending at F. These results rely on the assumption w < 1, see (JT]). 



The Bianchi type II subset 

Let us consider the Bianchi type II subset Bn^ given by N a — Np = 0, 7V 7 > 0. In this 
case, the 7-direction is singled out, while there exists a discrete symmetry associated with the 
interchange of the a- and /3-direction. The LRS subset E Q = E^ is a subset of codimension one 
which divides the state space into two equivalent parts, the subsets {E Q > E^} and {E Q < E^} 
(related by permuting the a- and /3-axes). 
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Since the past singularity is of particular interest in our considerations, it is convenient to 
choose the time direction towards the past, i.e., to use a reversed time variable r_ according to 

T_ = — T (18) 

which we do in the remainder of this section; accordingly, approach to the past singularity 
means t_ — * oo. (In Section [5] we will see that vacuum type II orbits are the building blocks 
for the Mixmaster map and the closely related Kasner map (BKL map) . It is a well established 
convention that forward iterations of these maps are directed towards the singularity. To agree 
with this convention the use of a past-directed time variable in a discussion of Bianchi type II 
models thus suggests itself.) 

On Bn 1 , the Gauss constraint E 2 + i JV2 +0 = 1 can be used to replace N 7 by as a dependent 
variable. The system ([7]) thus becomes 

^ = (2 - ? )E a//J + 3 S a/f3 , ^ = (2 - q)Z y +%, ^- = -SI [2q (1 + 3w)] , (19) 

where q = 2£ 2 + |(1 + 3w)Sl and 3 S a/p = -4(1 - E 2 - Si), % = 8(1 - E 2 - SI); we have 
£ 2 + 0<l. 

Let us first consider the vacuum subset B]f c ', i.e., = 0. There do not exist any fixed points 
in the type II vacuum subset BJf c ', but the boundary of the vacuum subset coincides with 
the Kasner circle K°. The orbits of (fT9|) form a family of straight lines in B™ c ', where each 
orbit connects one fixed point on K° with another fixed point on K°; hence each orbit is 
heteroclinic [ID] . Following the nomenclature of [30] we call these heteroclinic orbits Bianchi 
type II transitions, because each orbit can be viewed as representing a transition from one 
Kasner state to another. We denote these transitions by T/v^ , where each Tm 1 emanates from 
{"fafi) U Q 7 U (j(3a). If the initial point is a point of sector (7a/?), then the final point is a point 
of (cry/3) U {Q Q } U (a[3~f); interchanging a and j3 yields the transitions emanating from (p/(3a); 
if the initial point is Q 7 , the final point is T 7 ; the points T Q and Tp are not connected with 
any other fixed point (they are 'fixed points' under the present 'type II map'), see Figure [3] 

Let (SJj, £jg, E 7 ) = (3pjj — 1, 3p^— 1, 3p 7 — 1) denote the initial fixed point on K°; this point can 
be represented in terms of the Kasner parameter u — u 1 by using ([17]) . The orbit (transition) 
emanating from this fixed point is given in terms of an auxiliary function r/ = t)(t-) by 

E a//3 = 2[l-7 ? ]+7 7 S 1 Q//3 , E 7 = -4[l-7 ? ]+7 ? E 7 , (20a) 

where 77 is determined by the equation 

■p- = 2(1 - E 2 )r, with (1-E 2 ) = ~(9-v)(v-l) and g = \ + - ± U \ (20b) 

and the conditions that lim T _^_oo rj = 1 and lim T _^ +oc rj = g. The quantity g is in the interval 
(1,3) for u G (l,oo); u = 1 corresponds to g = 3 and describes the orbit Q 7 — > T 7 ; u = 00 
corresponds to g = 1 and describes the 'isolated' points T Q , Tg. Since rj increases from 1 to 
g as r_ goes from —00 to +00, g is called the growth factor [3D]. In Section [5] the transitions 
Tn , as represented by (|2"Pj) . will appear as the building blocks for the Mixmaster/Kasner map. 

While there do not exist any fixed points in the vacuum subset of S/v^, there exists one fixed 
point in £>jv t with Q > 0, the Collins-Stewart fixed point CS 7 , which corresponds to one repre- 
sentation of the LRS solutions found by Collins and Stewart [38]. CS 7 is given by (£ a , E^, E 7 ) = 
|(1 + 3w)(l,l,-2) and Cl = 1- i(l + 3w) (which yields 7V 7 = | yT^w^fl + 3w ). The fixed 
point CS 7 is the source (w.r.t. r_) for all orbits in Bn with Q > 0. In the limit t_ — > 00 all 
solutions in £>at 7 \{CS 7 } converge to fixed points on the Bianchi type I boundary of 2?/v T : There 
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Figure 3: The type II transitions 7/Vi on the Bn 1 subset; by definition, Ni ^ along T^, 
while N2 — N3 — 0. The projections of these transition onto (Si, E2, E3)-space are straight 
lines, which possess a common focal point Mi characterized by (Ei,E2,E3) = (—4,2,2). The 
transitions T/v 2 , Tjv 3 on the subsets Bn 2 , Bn 3 are obtained by permutations of the axes, see 
Figure [SI The arrows indicate the direction of time towards the past. 

exists one orbit (which corresponds to an LRS solution) that converges to F as t_ — > 00; every 
other orbit converges to a fixed point on (aj(3) U {Q a } U (a(3-f) or {(3"fa) U {Q/3} U (/3a-f) on 
K°, or to T 7 (in the LRS case). For a detailed discussion of these results see [TO] . 

The Bianchi type VII subset 

In anticipation of Theorem 16.11 which implies that generic orbits of Bianchi type IX do not 
have a-limi^ points (w.r.t. the standard future directed time variable r) on any of the Bianchi 
type VIIo subsets Bn 1 n 2 , Bn 2 n 3 , Bn 3 N!, we refrain from giving a detailed discussion of these 
subspaces here. (However, note that in order to prove Theorem 16.11 a detailed understanding 
of solutions of Bianchi type VIIo is essential; in fact, in the proof of Theorem 16.11 Bianchi 
type VIIo is ubiquitous; we refer to [26] and [27].) In the present context it suffices to note that 
on each subset Bn^n? there exists a line of fixed points TL 7 given by (E Q , E^, E 7 ) = (—1,-1, 2) 
and N a = Np (so that f2 = 0). Since TL 7 emanates from the point T 7 S K° we call it the 
'Taub line.' Like T 7 itself, each of fixed points on TL 7 is associated with a representation of 
Minkowski spacetime (in an LRS type VIIo symmetry foliation). 

4 Asymptotic self-similarity 

In the previous section we have given the fixed points associated with the system ([7]) on £?ix- A 
local dynamical systems analysis of the fixed points shows whether or not these points attract 
type IX orbits in the limit r — > — 00 We merely state the results here and refer to [27] for 
details. 

5 For a dynamical system on a state space X, the o-limit set a(x) of a point x 6 X is defined as the set of all 
accumulation points towards the past (i.e., asr-t —00) of the orbit 7(1") through x. The simplest examples of 
oj-limit sets are fixed points and periodic orbits. 

"In this section we adapt to 10 26] and use the normal future-directed time variable r. 



4 ASYMPTOTIC SELF-SIMILARITY 



11 



K° Each fixed point K on K°\{Ti, T2, T3} is a transversally hyperbolic saddle (one stabk|l| 
mode, see Figure IU and three unstable modes; the stable manifold of K coincides with a 
vacuum type II transition orbit, see Figure[3]). The Taub points {Ti,T2,T3} are center 
saddles with a two-dimensional unstable manifold and a three-dimensional center mani- 
fold, where T a is excluded as an a- limit set on the center manifold, see |27j : consequently 
there do not exist any type IX solutions that converge to any of the points on K° as 

T — ► —OO. 

F The fixed point F on Bq is a hyperbolic saddle (with B*' as a two-dimensional stable 
manifold and an additional three-dimensional unstable manifold) . Accordingly, F attracts 
a two-parametric family of type IX orbits as r — * — 00. These solutions have a so-called 
isotropic singularity. 

CS a The fixed points CS a (a — 1,2,3) on B^ are hyperbolic saddles (with B^j as a three- 
dimensional stable manifold and an additional two-dimensional unstable manifold). The 
unstable modes are associated with the equations Np N'p \cs a = f (1 + 3io) (for [3 ^ a). 
Therefore, each of the fixed points CS Q attracts an (equivalent) one-parameter set of 
type IX orbits in the limit r — > —00. 

TL Q Each fixed point on TL„ on B]^' N ^ is a center saddle. On the three-dimensional center 
manifold (which coincides with BJJ^y , a ^ / 7 / a) the point acts as a (non- 
hyperbolic) sink, see |27j . Since there is a two-dimensional unstable manifold, there exists, 
for each fixed point on TL a , a one-parameter family of type IX orbits that converges to it 
as r — > —00; these orbits correspond to LRS solutions. (Conversely, generic LRS type IX 
solutions converge to TL Q , see, e.g. ,[10].) 




Figure 4: This figure depicts the Kasner circle K° and the stable variables for each sector; N a is 
the stable variable in sectors (a/37), ( a 7/3), and at the point Q Q . Expressed in the time variable 
r_, which is directed towards the past, these variables are the unstable modes. At a given fixed 
point, the associated unstable manifold orbit is a Bianchi type II transition, see Figure [3J 

Collecting the results we see that the solutions whose a-limit is one of the fixed points form a 
subfamily of measure zero of the (four-parameter) family of Bianchi type IX solutions. Following 



7 When we reverse the direction of time, i.e., when we use r_ instead of r, we must replace 'stable' by 
'unstable'. 
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the nomenclature of [26] we thus refer to these solutions as non-generic solutions of Bianchi 
type IX. Alternatively, to capture the asymptotic behavior of these solution, we use the term past 
asymptotically self-similar solutions. (Since a fixed point in the Hubble-normalized dynamical 
systems formulation corresponds to a self-similar solution, see e.g. [10] , solutions that converge 
to a fixed point are asymptotically self-similar.) 

The past asymptotically self-similar solutions comprise the LRS Bianchi type IX solutions. As 
seen above, generic LRS solutions converge to TL Q towards the past (and each solution that 
converges to TL a is LRS), but there exist exceptional LRS solutions that converge to F or CS a . 
The remaining orbits whose limit point is either F or CS Q correspond to past asymptotically 
self-similar solutions that are non-LRS. Clearly, every solution that converges to F or CS a is a 
non- vacuum solution, since f2 ^ at F and CS Q . 

It is natural to ask how the non-generic orbits are embedded in the state space -Bix- The LRS 
orbits form the three LRS subsets £RS a , which are the hyperplanes given by the conditions 
E/3 = S 7 , Np = 2V 7 , where (a/37) € {(123), (231), (312)}. The orbits whose a-limit set is the 
fixed point CS a (for some a) form the set CS a in Bix', we call CS a the Collins-Stewart manifold. 
The local analysis of the fixed point CS Q and the regularity of the dynamical system ([7]) imply 
that the Collins-Stewart manifold CS a is a two-dimensional surface; it can be viewed as a two- 
dimensional manifold with boundary embedded in -Bix (where this boundary corresponds to 
an orbit in -Bvii ) ■ Analogously, the orbits whose a-limit set is the fixed point F form the set 
1 in B\x, which we call the isotropic singularity manifold, since solutions converging to F are 
those with an isotropic singularity. The isotropic singularity manifold $ is a three-dimensional 
hypersurface; it can be viewed as a three-dimensional manifold with boundary. 

Generic Bianchi type IX models are those that are not asymptotic self-similar and thus con- 
stitute examples for asymptotic self-similarity breaking; for other such examples, see [3U1 HU] , 
The central theme in this paper is the past asymptotic behavior of the generic models. 

5 Facts about the Mixmaster and Kasner map 

The Mixmaster attractor Aix (alternatively referred to as the Bianchi type IX attractor) is 
defined to be the subset of Bix given by the union of the Bianchi type I and II vacuum 
subsets, i.e., A\x — -Bj ac ' U BJf c \ Since the type II vacuum subset consists of three equivalent 
representations we obtain 

Aix = K° U U B v ^ U BJ£- . (21) 

In this section we investigate the structures that the flow of dynamical system (J7|) induces on 
Aix- In particular we discuss the Mixmaster map, the Kasner map, and the era map. To agree 
with the well-established convention for these maps, the direction of time will be taken towards 
the past. 

The Mixmaster, Kasner, and era maps 

In Section [3] we have seen that the vacuum type II orbits, i.e., the orbits on S™ ', a = 1,2,3, 
are heteroclinic orbits that emanate from and converge to fixed points on the Kasner circle K° , 
see Figures [3] and [5] In accord with Section [3] we refer to these orbits as transitions and denote 
them by Tjv Q , a — 1, 2, 3. 

The type II transitions are the building blocks for the analysis of the Mixmaster attractor. 
By concatenating transitions we obtain a sequence of transitions, also known as a heteroclinic 
chain. Since each fixed point on K° (except for the Taub points) is the initial value for one 
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Figure 5: Projection of the type II transitions Tm 1 Tn 2 , T/v 3 on the type II subsets B™^, 
£>™ 2 C ', /8™ 3 C ' onto (Ei, E 2 , E 3 )-space. Concatenation of these transition orbits yields sequences 
of transitions. The arrows indicate the direction of time towards the past. 



single transition, sec Figures [4] and the concatenation of transitions is unique: Each fixed 
point (except for the Taub points) generates a unique sequence of transitions. Note, however, 
that the 'direction of time' is relevant. For each fixed point P on K°, which is not one of the 
Taub points, there exists one single transition emanating from P at t_ = — oo, but there are two 
transitions converging to P as t_ — > oo. Therefore, concatenating transitions in the reversed 
direction of time leads to ambiguities. (In terms of the standard future-directed time variable 
r we have the converse statement: It is possible to make unambiguous retrodictions, but not 
predictions.) 

Let I = 0, 1, 2, . . . and let P; G K° denote the initial point of the transition (P/ is also the 
end point of the (I — 1)^ transition). We refer to the sequence (P/)/ e N of Kasner fixed points, 
which is induced by the sequence of transitions, as being generated by the Mixmaster map. The 
Mixmaster map can be visualized by a map in (Ex, E 2 , E 3 )-space, obtained by inscribing K° 
in a triangle with corners at (E 1; S 2 , E 3 ) = (—4, 2, 2) and cyclic permutations, from which the 
(projections of the) transition orbits 'originate' as straight lines; see Figure 




Figure 6: Concatenating type II transition orbits we obtain sequences of transitions — 
heteroclinic chains. The discrete map governing the associated sequence of fixed points on 
K° is the Mixmaster map. The arrows indicate the direction of time towards the past. 

Let the initial Kasner state of a transition be represented by the Kasner parameter u = u\ 
where we assume u 1 < oo, since neither of the Taub points Ti, T 2 , T 3 can be the initial value 
for a transition. Inserting (|15[) and (fT7|) into (f2T)]) we find that a transition maps the Kasner 
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parameter u 1 to the parameter u f , where 

{ u 1 -1 if v} £ [2, oo) . 

« f = S • , . (22) 

\ K-l)" 1 if «'€ [1,2]. 

The information contained in this Kasner map suffices to represent the collection of all transition 
orbits (as a whole). However, for each particular mapping u 1 h- > u f , there exist six (equivalent) 
associated transitions jf| this is simply because u 1 characterizes the initial Kasner fixed point on 
K° only up to permutations of the axes. Hence, in order to reconstruct a particular transition 
from (|22p , we supplement (|22|) with information about the initial sector of the transition, which 
determines the position of the axes. 



In terms of the Kasner parameter, a sequence of transitions corresponds to an iteration of {22]). 
Let Z = 0,1,2,.. . and let ui denote the initial Kasner state of the transition. This transition 
maps ui to it;+i, i.e., 

Jtb transition f U l ~ 1 if U l G [ 2 , Oo) , 

Mi * 1 Ul + 1 = 1 / ,x_1 .p r, «vi ( 23 ) 




if uiG[l,2] 

We refer to this map as the (iterated) Kasner map (which is also known as the BKL map [5]). 
Since each value of the Kasner parameter u € (l,oo) represents an equivalence class of six 
Kasner fixed points, the Kasner map can be regarded as the map induced by the Mixmaster 
map on these equivalence classes via the equivalence relation. 

In a sequence (ui){=o,i,2,... that is generated by the Kasner map (|23p . each Kasner state it; 
is called an epoch. Every sequence (it;)z=o,l,2,... possess a natural partition into pieces (which 
contain a finite number of epochs each) where the Kasner parameter is monotonically decreasing 
according to the simple rule ui i— ► ui+\ = ui — 1; these pieces are called eras [2]. An era begins 
with a maximal value ii; in of the Kasner parameter (where ui ia is generated from wz in -i by ui in = 
[w; in _i — 1] ), continues with a sequence of Kasner parameters obtained via ui t— > ui + \ = — 1, 
and ends with a minimal value ui mt that satisfies 1 < it; out < 2, so that u; out +i = [u; out — 1] _1 
begins a new era. 

6.29 -> 5.29 -> 4.29 -> 3.29 2.29 -> 1.29 -> 3.45 -> 2.45 -» 1.45 -> 2.23 -> 1.23 -» 4.33 

era era era era 

(24) 

Let us denote the initial (= maximal) value of the Kasner parameter u in era number s (where 
s = 0, 1, 2, . . .) by u s . Following [7] we decompose u s into its integer part k s — [u s ] and its 
fractional part x s = {u s }, i.e., 

u s = k s + x s , where k s = [u s ] , a; s = {u s } . (25) 

The number k s represents the (discrete) length of era s, which is simply the number of Kasner 
epochs it contains. The final (= minimal) value of the Kasner parameter in era s is given by 
1 + x s , which implies that era number (s + 1) begins with 

1 _ 1 

u s+l — — 7 T" • 

The map u s i— > u s+ i is (a variant of) the so-called era map; starting from uo = uq it recursively 
determines u s , s — 0, 1, 2, . . ., and thereby the complete Kasner sequence (wi)i=o,i,...- 



8 In the exceptional case u 1 = 1 there exist only three (equivalent) associated transitions: Qi — > Ti, Q2 — > T2, 
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The era map admits a straightforward interpretation in terms of continued fractions. Consider 
the continued fraction representation of the initial value, i.e., 

u = k H = [fc ; ki, k 2 , fc 3 , . . .] . (26) 

h 



k 2 + --- 

The fractional part of uo is xq = [0; k\, k 2 , k 3l . . .]; since Ui is the reciprocal of xq we find 

u 1 = [k 1 ;k 2 ,k 3 ,k A ,...}. (27) 
Therefore, the era map is simply a shift to the left in the continued fraction expansion, 

u s = [k s ;k s+1 ,k s+2 , ■ ■ ■} >-> u s+1 = [k s+1 ;k s+2 ,k s+3 , . . .} . (28) 

The properties of the Kasner sequence depend on the initial value Uo = uq. 

(i) If and only if the initial Kasner parameter uq is a rational number, i.e., if and only if 
uo S Q, then its continued fraction representation is finite, i.e., 

u = [ko;ki,k2,...,k n ] , (US) 

where k n > 1. Therefore, there exists only a finite number of eras (where the last one 
begins with u„ = k n ), and the Kasner sequence is finite. At the end of era number n, 
the Kasner parameter reaches u — 1, which subsequently terminates the recursion (|23p . 
Since Q is a set of measure zero in R, this case is non-generic. 

(ii) A quadratic irrational (quadratic surd) is an algebraic number of degree 2, i.e., an irra- 
tional solution of a quadratic equation with integer coefficients. If and only if the initial 
Kasner parameter u is a quadratic irrational, i.e., if and only if uq = q\ + ^fq-2., where 
qi € Q and q 2 € Q is not a perfect square (i.e., y/q^ £ Q), then its continued fraction 
representation is periodic, i.e., 

u = [fc ;fci,...,An, ; (EHU) 

the notation is such that the part in parenthesis, i.e., (fci, . . . , kn), is repeated ad infini- 
tumj^l Consequently, the era map becomes periodic (after the rfi 1 era), and we thus obtain 
a periodic sequence of eras and a periodic Kasner sequence (wz)z=o,i.2,...- It is straight- 
forward to see that while the period of the era sequence is n, the period of the Kasner 
sequence is (k\ + • • • + k n ); see the examples below. Since the set of algebraic numbers of 
degree two (or cquivalently the set of equations with integer coefficients — it is a subset of 
N 3 ) is a countable set, case (ii) is also non-generic. 

(iii) An irrational number is called badly approximable if its Markov constant^ is finite. If and 
only if uo is badly approximable, then the coefficients (partial quotients) in its continued 
fraction representation are bounded, i.e., 

uo = [ko; k\, ki, fe, • • .] with ki < K V« (|29lii) 

for some positive constant K . Consequently, the sequence of eras and the Kasner sequence 
(ui)/=o,i,2,... are bounded, i.e., ui < K V7. Obviously, case (ii) is a subcase of case (iii). 



9 If uo = 9i + y/q2 > 1 is a quadratic irrational such that qi — v / <?2 £ (—1, 0), then uo = Uki, . . . , kn)] , i.e., 
the continued fraction is purely periodic without any preperiod. 

10 For i£8, let denote the distance from x to the nearest integer, i.e., = min n gz |x — n|. The Markov 
constant M(x) of a number x £ M\Q is denned as M(x)~ 1 = liminfj} 3lwoo n ||nx||, see 141 1 . It is known that 
M(x) > x/5 for all x e 
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(Note, however, that there is probably no relationship between case (iii) and algebraic 
numbers of degree greater than 2, since it is expected that these numbers are well ap- 
proximate.) The set of badly approximate numbers are a set of Lebesgue measure zero, 
hence this case is non-generic. 

(iv) If and only if the initial Kasner parameter uo is a well approximable irrational number, 
then the partial quotients ki in the continued fraction representation 



uo = [k Q ; ki, k 2 , k 3 , . . .] 

are unbounded (and we can construct a diverging subsequence from the sequence of partial 
quotients (fci)igN). This is the generic case, and hence generically the Kasner sequence 
(ui)i—o t x t 2 t ... is infinite and unbounded. 

In terms of continued fractions, the Kasner sequence (wi)z e R generated by uq = [fcoj k\, fc 2 , . . . ] 
is 

u = u = [fc ; fci,fc 2 , ...] -> [k - 1; fci, k 2 , . ■ .] -> [k - 2; fci, fc 2 , . . .] ->...—> [l; fci, fc 2 , . . .] 
— > ui = [fci; fc 2 , fc 3 , . ..] — > [fci - 1; fc 2 , fc 3 , . . .] — > [fci - 2; fc 2 , fe, ■ ■ • ] —►...—> [l; fc 2 , &3, . . . ] 
— * u 2 = [fc 2 ; fc 3 , fc 4 , . . . ] — * [fc 2 - 1; fc 3 , fc 4 , . . . ] — * [fc 2 - 2; fc 3 , fc 4 , . . . ] — * . . . 

Let us give some examples for periodic era sequences and Kasner sequences. If Uq = [(1)] 
; I + y/S)/2, which is the golden ratio, then u s = (1 + \/5)/2 Vs. It follows that the Kasner 
sequence is also a sequence with period 1, 

(ui)i eN : |(l + \/5) -> |(1 + V5) -» i(l + V5) ^ i(l + V5) ^ i(l + V5) ^ ... 

If w = [(2)] = 1 + V2, then u s = 1 + y/2 Vs; hence the era sequence is a sequence of period 1. 
However, the associated Kasner sequence has period 2, 



(1 + V2) -► V2 -h. (1 + V2) ^ V2 -h. (1 + V2) -* V2 -* (1 + V2) -* V% ■ 



Analogously, the initial value uo = [(3)] = (3 + vT3)/2 generates an era sequence of period 1 and 
an associated Kasner sequence of period 3. Finally let uq = [(2, 4)] = 1 + v/3/2 ~ 2.2247. Then 
the era sequence has period 2 with u„ = [(2, 4)] ~ 2.2247 for even n and u„ = [(4, 2)] ~ 4.4495 
for odd n. The associated Kasner sequence (u;)/gN has period 6, 

2.22 -> 1.22 -> 4.45 -> 3.45 -> 2.45 -> 1.45 -> 2.22 1.22 -> 4.45 -> 3.45 2.44 -> 1.45 



In the state space description of sequences (in terms of the Mixmaster map), an epoch is simply 
a point P; on the Kasner circle. (It is one of the six points in the equivalence class associated 
with the Kasner parameter «;.) Transitions connect epochs and thus generate the Mixmaster 
map. 

The Kasner parameter u can be employed to measure the (angular) distance of a point P on 
K° from the Taub points or from the non-flat LRS points: If u ^> 1, then P is at an angular 
distance of approximately it" 1 from one of the Taub points. On the other hand, in the vicinity 
of the non-flat LRS points (where u— 1<1),m- lisa linear measure for the angular distance 
of P from the closest of the non-flat LRS points. If u < 2, then P is closer to one of the non-flat 
LRS points Q a than to any of the Taub points T a . Therefore, in the state space picture, an 
era can be described as a (finite) sequence of points K° 3 P;, l- m < I < Z out , obtained from 
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the Mixmaster map, whose angular distance from the Taub points is monotonically increasing. 
The era begins with a fixed point P; in that is close to one of the Taub points T a (the preceding 
point P/ in -i was closer to one of the other two Taub points T^, /3 ^ a, than to Tq,); then 
the distance from the Taub point T a monotonically increases until, for P; out , it exceeds the 
distance to the non-flat LRS points; this is the terminal point for the era; it is connected by 
the following transition with the initial point of the next era; a good illustration is Figure [7(d) | 
where each era consists of three epochs. Note that due to the equivalence of Kasner points, 
there exist several realizations of one and the same Kasner sequence as Mixmaster sequences 
in the state space picture; this is exemplified by the heteroclinic cycles in Figure [7(a)1 and [7(b)] 

The state space description of the cases (i)-(iv), which are characterized by the initial Kasner 
parameter uq — uq, is the following: 



(i) Iff uo G Q, see (|29ip , the Mixmaster sequence of Kasner fixed points (P;);=o.i,... is finite. 
After a finite number of transitions, at the end of era n, the sequence reaches one of the 
LRS points Q a (where u = 1); a last transition follows, namely the transition Q Q — > T Q , 
and the sequence terminates in one of the Taub points. 

(ii) Iff uo is a quadratic surd, i.e., uq = q\ + */#2 for some q\, q% G Q where q<i is not a perfect 
square, see (|29ti[) . then the Mixmaster sequence of Kasner points (P;); £ n is eventually 
periodic, where the period is a multiple of (ki + ■■■ + kn); see Figure Viewed as a 
periodic sequence of transitions (which are heteroclinic orbits) we obtain a heteroclinic 
cycle. In Figure[7]we give some of the heteroclinic cycles associated with Kasner sequences 
with periods 1, 2, and 3 in their projection onto (Si, E2, £3)-space. Note that due to 
permutation symmetry there are several cycles associated with a given periodic Kasner 
sequence. 

(iii) Iff uo is a badly approximable irrational number, see (|29iiip . there exists a neighborhood 
of the Taub points T Q such that the Mixmaster sequence of Kasner points (P;)/ e N does 
not enter this neighborhood. This is simply because there exists a maximal value of the 
sequence (u ; ) ;eN . 

(iv) Iff uo is a well approximable irrational number, see (|29ivp . then the Mixmaster sequence 
(P/)i S N comes arbitrarily close to the Taub points. This is the generic case. 



In the following we analyze the generic case (iv) in more detail. Let uq be a well approximable 
irrational number, i.e., a number whose continued fraction expansion 

u = [k ;ki,k 2 ,k 3 , . . .] (30) 

defines an unbounded sequence (fci)ieN- By construction, era number i contains ki epochs 
(which we call its length). A natural question to ask concerns the distribution of the partial 
quotients ki. For a 'typical' well approximable irrational number, how often does the number 
1 appear in the sequence (ki)i^iq? How often the number 2? And what about the number 
1000? The answer is given by Khinchin's law [32] • Let P n (k=m) denote the probability that 
a randomly chosen partial quotient among . . . , k n ) equals m G N. In the asymptotic limit, 
i.e., for P(k=m) — linin^oo P n (k=m) we have 

P{k =m) = 2 log p±±) - 2 log (_^_) , (31) 
v m + 2 m + 1 

i.e., the partial quotients of the continued fraction representation of (|30p are distributed like 
a random variable whose probability distribution is given by (|3ip . (By (|3ip . the number 1 
appears in 42% of the slots, the number 2 in 17%, and the number 1000 in 1.4 * 10~ 6 % of the 
slots.) Khinchin's law applies for almost all numbers Uq. 
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(a) One of the two heteroclinic cycles 
associated with the Kasner sequence 
generated by uq = [(!)]• 




(c) One of the heteroclinic cycles asso- 
ciated with the Kasner sequence gener- 
ated by uq = [(2)]. 




(b) The second of the two heteroclinic 
cycles associated with the Kasner se- 
quence generated by uq = [(!)]• 




(d) One of the heteroclinic cycles asso- 
ciated with the Kasner sequence gener- 
ated by uq = [(3)]. 



Figure 7: Examples of heteroclinic cycles associated with era sequences with period 1. The 
Kasner sequences have period 1, 2, and 3, respectively; the period of the heteroclinic cycles is 
a multiple of that period. Note that the direction of time is towards the past. 



For a (generic) Kasner sequence (u;)/ 6 n with initial parameter uq = [ko; hi, fe, • ■ •] and its 
associated era sequence (u s ) se rsj, where u s = [k s ; fc s +i, fc s +2, • ■ ■], the expression P{k=m) of (|3Tj) 
represents the probability that a randomly chosen era of (u s ) s£ n has length to; this corresponds 
to the probability that the initial value of an era is contained in the interval [to, to + 1). In this 
manner, the probability distribution (|31[) makes possible a stochastic interpretation of generic 
Kasner sequences. 

The probability distribution (13ip results in extraordinary properties of the Mixmaster/Kasner 
map, which will be of crucial importance in the considerations of Section [71 We do not discuss 
details here but refer to future work; however, we cannot refrain from giving a teaser: For a 
generic Kasner sequence (ui)jeN and its associated era sequence (u s ) se N there exist infinitely 
many eras such that the length (i.e., the number of epochs) of the era is larger than nlogn; 
however, for sufficiently large n, the length is guaranteed to be bounded by nlog 2 n. (For a 
proof of this result by Borel and Bernstein see [43]; see also [44].) Properties of this kind 
underline the remarkable intricacies of the heteroclinic structures on the Mixmaster attractor. 
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6 Mixmaster facts 

In this section we turn to what is known about the past asymptotic dynamics of generic type IX 
solutions. The main Mixmaster fact is Ringstrom's 'Bianchi type IX attractor theorem'. 

The main Mixmaster fact 

Consider a solution of Bianchi type IX that is cither vacuum or associated with a perfect 
fluid satisfying —^<w<l. Recall that such a solution is called generic if it is not past 
asymptotically self-similar, i.e., if its a-limit set is neither the point F, nor any of the points 
CS a , nor a point on TL Q ; in other words, a generic solution corresponds to an orbit in Bix 
that is neither contained in 5F, nor in CS a , nor in CRS a . Therefore, the set of generic Bianchi 
type IX states is an open set in £?ix- (To conform with [TUM^S] we use the future directed time 
variable r.) 

The main results concerning generic Bianchi type IX models are due to Ringstrom [26]; these 
results rest on earlier work that is reviewed and derived in [10] , and on [24j, [25] . In the following 
we state the main theorem in a version adapted to our purposes. 

Theorem 6.1 ([IS]). Let (Si, S2, S3, N±, N2, Ns)(t) be a generic solution of Bianchi type IX, 
i.e., a generic solution of (J7J) in Bix- Then 

A n = N X N 2 + N 2 N 3 + N 3 N t -> and fl -> (32) 

as t — > —00. 

Note that this Theorem applies to both the fluid and the vacuum case; in the latter case ([32|) 
becomes An — > and il = 0. 

The proof of Theorem 16.11 given in [26] is delicate. In the first part it is proved that the a-limit 
set of a generic Bianchi type IX solution is non-empty and must contain a point on the Kasner 
circle. The main part of the proof deals with the fact that the function An(r) is in general 
not monotone. There exist times where An(r) increases (as r — > —00); the associated growth 
must therefore be controlled and shown to be negligible compared to the overall decrease in 
An. This is done by a careful analysis of the equations and (approximate) solutions. In [27] 
we give an alternative and relatively short and succinct proof of Theorem 16.11 which is based 
on an in-depth understanding of the hierarchical structure of the dynamical system j7]) (as 
represented by Figure !!). It is important to note, however, that either of the proofs fail in 
the other Bianchi types that are conjectured to exhibit an oscillatory approach towards the 
singularity, i.e., the proofs fail for types VI_i/g and VIII. This is unfortunate, since there are 
reasons to believe that these models are more relevant than type IX as regards the dynamics 
of generic (inhomogeneous) cosmologies, see |27) . 

Using the concept of the Mixmaster attractor, cf. ([2"Tj) , we obtain an equivalent formulation of 
Theorem 16. II Let X(r) = (Si, S2, S3, Ni, N2, N 3 )(t) be a generic solution of Bianchi type IX. 
Then 

\\X(t) - Avt\\-> (r--oo), (33) 
where the distance \\X — -4ix|| is given as miny g _4 IX \\X — Y\\. 

Theorem 16.11 thus states that the attractor of generic type IX solutions resides on Aix', how- 
ever, whether the past attractor is in fact _4ix or merely a subset thereof remains open. (The 
terminology 'Mixmaster attractor' is seductive but might turn out to be quite misleading.) Like- 
wise, the theorem does not provide any direct information about the details of the asymptotic 
behavior of solutions. 
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Consequences 

The prerequisite for a deeper understanding of the asymptotic behavior and the oscillatory 
dynamics of generic type IX solutions is an understanding of the Mixmaster attractor. In 
Section [5] we have identified the structures on the Mixmaster attractor Aix that are induced by 
the flow of the dynamical system: heteroclinic cycles [case (ii)] and finite [case (i)] and infinite 
heteroclinic sequences [cases (iii) and (iv)]. All these structures qualify (a priori) as possible 
a-limit sets of generic type IX orbits. In conjunction with these results, the main theorem 16.11 
implies a number of further facts about the attractor — 'Mixmaster attractor facts', which we 
give as a list of corollaries. (For proofs see [27], and also [26].) 

1. A generic, i.e. not past asymptotically self-similar, type IX orbit possesses an a- limit 
point on Aix- 

2. If P e Aix is an a-limit point of a type IX orbit, then the entire heteroclinic cycle/sequence 
(Mixmaster sequence) through P must be contained in the a-limit set. 

3. If one of the Taub points {Ti,T2,T3} is an a-limit point of a type IX orbit, then the 
a-limit set contains Kasner fixed points associated with arbitrarily large values of the 
Kasner parameter u. 

4. For generic solutions of Bianchi type IX the Weyl curvature scalar C a bcdC ahcd (and there- 
fore also the Kretschmann scalar) becomes unbounded towards the past. 

5. Taking into account both the expanding and contracting phases of Bianchi type IX so- 
lutions, generic Bianchi type IX initial data generate an inextendible maximally globally 
hyperbolic development associated with past and future singularities where the curvature 
becomes unbounded. 

6. Convergence to the Mixmaster attractor is uniform on compact sets of generic initial data: 
Let X be a compact set in _E?ix that does not intersect any of the manifolds J, CS a , £RS a , 
so that each initial data x £ X generates a generic type IX solution. Let X(x; r) denote 
the type IX solution with X(x, 0) = x. Then 

\\X(x;t)-Aix\\^0 (r-f-oo) (34) 

uniformly in x £ X. 

Corollary 3 implies that the a-limit set contains an infinite set of Kasner fixed points in a 
neighborhood of the Taub point (s), but this set is not necessarily a continuum of fixed points; 
cf. the previous discussion about possible a-limit sets. If a Kasner point with u £ Q is contained 
in the a-limit of an orbit, so is a Taub point. This is an immediate consequences of the results 
of Section [5j case (i) . On the other hand, if the a-limit set of an orbit is a heteroclinic cycle 
or a heteroclinic sequence (possibly in combination with a cycle) associated with cases (ii) and 
(iii) of Section [5j then there exists a neighborhood of the Taub points whose intersection with 
the a-limit set is empty. 

Note that the oscillatory behavior of asymptotic type IX dynamics, which we unfortunately 
know no details about, constitutes an example of asymptotic self-similarity breaking [39] . In 
order to make progress as regards the details of the asymptotic oscillatory behavior, it is natural 
to first establish the connection between the Mixmaster/Kasner map and dynamics for a finite 
time interval At, discussed next. 
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Finite Mixmaster shadowing 

To make contact between the Mixmaster map and Bianchi type IX asymptotic dynamics we 
introduce the concept of finite Mixmaster shadowing which formalizes the following basic idea: 
Given a sequence of transitions we can choose type IX initial data sufficiently close to the initial 
data of the sequence so that the type IX solution generated by this data remains close to the 
sequence for some 'time'. 

Let Po be a Kasner fixed point (but Po ^ {Ti, T2, T3}) and let uq be the associated value of 
the Kasner parameter. There exists a unique sequence of transitions (7i)i£fs and an associated 
Mixmaster sequence (P;)z e N with Po as initial data; the associated Kasner sequence is (ui)i^. 
(If ito G Q, the sequence terminates at one of the Taub points after a finite number of transi- 
tions.) We shall make the definition that a type IX solution shadows a finite piece (^)/=o,i,...,i 
of the sequence of transitions if it is contained in a prescribed (small) tubular neighborhood 
of (7~i)i=o,i,...,l- However, a standard e-neighborhood of the sequence fails to be a reasonable 
measure of closeness, because in the vicinity of a Taub point the transitions lie so dense that 
consecutive transitions are not separated from each other by their respective e-neighborhoods. 
Therefore, the introduction of adapted tubular neighborhoods is necessary to take into account 
the sensitivity of the flow at the Taub points and to capture more accurately the intuitive idea 
of shadowing. 

Let e > be small. A ' Taub-adapted neighborhood' of the Mixmaster sequence (P;); e pj is a 
sequence (]Ai)i<zn of open balls, where Ui is centered at the point P; and has radius eit ; -2 , i.e., 
Ui = {X G Bix ■ \\X -Pi|| < eu^ 2 }. (The radius eu l 2 is chosen to ensure that the intersection 
of the ball Ui with the Kasner circle induces more or less a standard e-neighborhood (it; — e, Ui+e) 
of the Kasner parameter ui; recall from Section [5] that u^ 1 measures the angular distance of 
a fixed point from a Taub point.) A Taub-adapted tubular neighborhood of the sequence of 
transitions (Ti)i & ^ is the sequence of tubes (Vz); 6 n that linearly interpolate between Ui and Ui + \. 
Based on this definition we say that a type IX solution X(t) shadows a finite piece ( / ?z);=o,i,...,l 
of the sequence of transitions if it moves in a prescribed Taub-adapted tubular neighborhood 
{7~l)l=o,i, i-e., if there exists a sequence of times {111)1=0,1 l+i such that X{r) is contained 
in Vi for all r G {ti+i, i~i\ for all < I < L. 

Making use of these concepts, a formulation of finite Mixmaster shadowing is the following: 
Let e > and LeR Consider the sequence of transitions {Ti)ien emanating from an initial 
Kasner point Po and its Taub-adapted tubular neighborhood (associated with e). Then there 
exists 5 e > such that each type IX orbit X{t) that is generated by initial data Xo with 
l|Ao — Poll < $e shadows the finite piece {Ti)i=o,i,...,l 01 the sequence of transitions. (A proof 
of this statement — in a slightly different form — has been given by Rendall [23] . Alternatively, 
one can invoke the regularity of the dynamical system, the center manifold reduction theorem 
and continuous dependence on initial data.) Evidently, 5 e depends on the choice of e and L. 
More importantly, however, S e depends on the position of Po — shadowing is not uniform; in 
particular, if we consider a series of initial points that approach one of the Taub points, then 
5 e necessarily converges to zero along this series — shadowing is more delicate in the vicinity of 
a Taub point. We will return to this issue in some detail in the next Section. 

Finite Mixmaster shadowing concerns any generic type IX orbit X{t). Let P G K° be an a-limit 
point of the type IX orbit X{t); without loss of generality we may assume that P is not one 
of the Taub points. (The existence of such a point is ensured by Corollaries HH3] of Section [6]) 
For simplicity we assume that P is associated with an irrational value of the Kasner parameter, 
which guarantees that the sequence (P;);epj emanating from P = Po is an infinite sequence. 
Since P is an a-limit point of X{t), there exists a sequence of times {r n ) n ^, r„ — ► —00 
(n — > 00), such that X{r n ) — > P {n — » 00). Therefore, we observe a recurrence of phases, where 
the orbit X{t) shadows (P;)i S N with an increasing degree of accuracy, i.e., shadowing takes 
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place for increasingly longer pieces of the sequence or in ever smaller neighborhoods. (If the 
Kasner parameter of P is rational, then the sequence (P/); 6 n is finite and terminates at a Taub 
point. X(t) will shadow this finite sequence recurrently with an increasing degree of precision.) 
We will return to this issue in some more detail in the subsection 'Stochastic Mixmaster beliefs' 
of Section [7] 

The concept of shadowing leads directly to the concept of approximate sequences which we 
introduce next. Consider a generic type IX orbit X(t) = (Si, £2, £3, Ni, N2, Ns)(t) and the 
function ||-X"(t) — K°||, where the distance \\X — K°|| is given as miny gK o \\X — Y\\. When the 
orbit X(t) traverses a (sufficiently small) neighborhood of a fixed point on K°\{Ti, T2, T3}, 
the function ||A(r) — K°|j exhibits a unique local minimum. This is immediate from the 
transversal hyperbolic saddle structure of the fixed point. (However, the flow in the vicinity 
of the Taub points is more intricate, since these points are not transversally hyperbolic.) It 
follows that the function \\X(t) — K°|| can be used to partition X(t) into a sequence of segments 
in a straightforward manner I 11 ! The local minima of ||X(t) — K°|| form an infinite sequence 
(T"i)ieN such that 77 — > — 00 as I — > 00. (This follows directly from Corollaries Q] and [2] because 
X(t) has a-limit point (s) on the Kasner circle and a- limit points on the type II subset.) A 
segment of X(t) is defined to be the solution curve between two consecutive minima, i.e., the 
image of the interval (77+1,77]. (Note that r — > —00 in the approach to the singularity, while 
the discrete 'time' I is past-directed, i.e., I — > 00 towards the singularity. This convention is 
chosen to agree with the standard convention for the Mixmaster and the Kasner map.) In the 
asymptotic regime, i.e., in the approach to the Mixmaster attractor, finite shadowing entails 
that a finite sequence of segments will resemble a finite sequence of type II transitions; this 
assumes, however, that the type IX orbit does not come too close to any of the Taub points. 
(In the neighborhood of a Taub point the flow of the dynamical system is much more intricate. 
This might yield recurring interruptions of the 'standard behavior'.) We call the type IX orbit 
in its segmented form an approximate sequence of transitions. 

In addition, we define a sequence of 'check points' that is associated with an approximate 
sequence of transitions. Each minimum 77 of ||^(t) — K°|| is associated with a Kasner fixed 
point P; (a 'check point') that is defined as the minimizer on K° of the distance between X(t{) 
and K°. (Note that the check points (P;);eN do not lie on the type IX orbit X(t).) Since 
each check point P; is associated with a value ui of the Kasner parameter, the sequence (P;); e N 
induces a sequence (m;); s n- The value ui + \ is in general not generated from ui by the exact 
Kasner map (|23[) , but differs from that value by an error of Sui . In our terminology, the sequence 
(P;); e N is an approximate Mixmaster sequence; (ui)i e m is an approximate Kasner sequence, see 
Figured 

The approximate Kasner sequence ui associated with a type IX orbit X(t) does not follow 
the Kasner map (|23p exactly. A natural question to ask, however, is whether the errors Siii 
of the approximate Kasner sequence converge to zero as I — > 00 or not. If Sui — > as / — > 00 
for a type IX orbit, this means that its dynamics is completely described by an 'asymptotic 
Mixmaster /Kasner map', i.e., by a map that converges to the Mixmaster/Kasner map towards 
the singularity. However, if 5ui -/* as / — * 00, then the evolution is interrupted repeatedly — 
infinitely many times — by phases where the dynamics is completely different from the Mixmas- 
ter dynamics, e.g., 'eras of small oscillations' [2]. (In the present state space description an era 
of small oscillations is associated with type VIIo behavior in the vicinity of one of the Taub lines 
TL Q , where Np and N 1 are small and of the same order; for details we refer to the discussion 
of type VII dynamics in [27].) In the subsection 'Stochastic Mixmaster beliefs' of Section[7]we 

llr The definition of a partition of X(t) into segments seems natural; it is important to note, however, that any 
definition depends on the formulation of the problem and is to a certain extent arbitrary. For instance, instead of 
using the minima of || X (r) — K^ 1 1| one might prefer to analyze the projection of the orbit onto (Si, £3, Ss)-space 
and use the extrema of S 2 (r). However, the conclusions drawn from any construction of segments are quite 
insensitive to the details of the definition. 
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Figure 8: A type IX orbit decomposes into segments X{n) — » X(ti + i). The points X{ji), 
I E N, are the local minima of the distance between X(t) and K°. The 'check point' P; is 
the point on K° that is closest to X{r) (for r in a neighborhood of 77). Since each check 
point P; is associated with a value iii of the Kasner parameter, the 'approximate Mixmaster 
sequence' (Pz);eN induces an 'approximate Kasner sequence' (ui)i^. In general, Pi + i/ui + i are 
not generated from Pi/ui by the exact Mixmaster/Kasner map (|23[) . but differ by an error 
SPi/Sii. 

will investigate the behavior of Sui along type IX orbits in detail. 

7 Mixmaster beliefs 
Attractor beliefs 

There remain several important open problems. In the following we will address the most 
pressing questions; in particular we will transform vague beliefs into refutable conjectures, and 
give arguments in their favor. 

An immediate question is the following: What are the actual a- limit sets of type IX orbits? 
Consider a type IX orbit 7. The a-limit set 0(7) contains a number of Kasner fixed points (and 
the associated type II transitions) , each of which is characterized by a particular value of the 
Kasner parameter u. The set of Kasner parameters obtained in this way from (2(7) we denote 
by C/(7). The question of which form 1/(7) can take for different orbits 7 is open, the most 
significant issues being the following: 

(i) Is it possible that there exists 7 whose limit set U (7) consists of rational numbers only? It 
is straightforward to exclude that 1/(7) coincides with Q itself (or a dense subset thereof), 
the reason being that a-limit sets are necessarily closed. Another a priori constraint is 
that U (7) must contain u — 00 (which characterizes the Taub points) if U (7) contains 
a rational number; see Section [Sj Corollary [3] of Section [5] then implies that ^(7) is 
unbounded, i.e., 1/(7) contains arbitrarily large values of the Kasner parameter. A set 
that is compatible with these basic requirements is, e.g., ^(7) = NU {00}. Whether there 
exist orbits 7 such that U (7) takes this (or a related) form remains open. 

(ii) Is it conceivable that there exist orbits 7 such that 0(7) is a heteroclinic cycle, see Fig- 
ure 0? In this case the set U{^) is generated by a quadratic surd u — {(k~i, fe, ■ ■ ■ ,k n )} 
via the Kasner map; see Section [S] in particular, f (7) is finite. However, whether orbits 
7 with this particular past asymptotic behavior really exist is an open problem. 
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(iii) Can there exist orbits 7 such that U fa) is bounded but contains infinitely many u-values? 
Is the Kasner sequence generated by a badly approximable number a candidate? The 
Kasner sequence (u/)z e N generated by a badly approximable number is an infinite sequence 
that is bounded. However, there must be at least one accumulation point of this sequence; 
if Ufa) contains (mz)/ 6 n, then Ufa) must contain the accumulation point of the sequence 
as well (since a-limit sets are closed). If this accumulation point is a well approximable 
number, Ufa) cannot be bounded; however, if the accumulation point is a quadratic surd, 
no inconsistencies arise. Hence, a priori, there might exist orbits 7 such that U fa) is not 
finite but still bounded. Whether this is indeed the case is doubtful, but hard to exclude 
a priori. 

An open problem that might be quite separate from the questions raised above concerns the 
behavior of all type IX orbits save a set of measure zero. 

Definition. The past attractor of a dynamical system given on a state space X is defined as 
the smallest closed invariant set A~ C X such that a{p) C A~ for all p 6 X apart from a set 
of measure zero J^5/ . 

Conjecture (Mixmaster attractor conjecture). The past attractor of the type IX dynamical 
system coincides with the Mixmaster attractor A\x (rather than being a subset thereof). 

Why is this a belief and not a fact? Theorem 16 . 1 1 implies that A~ C Aix', however, it is believed 
that A~ = Aix = K°UB V ^-UB V ^-UB V ^\ (The usage of the terminology 'Mixmaster attractor' 
for the set Aix reflects the strong belief in the Mixmaster conjecture.) It is difficult to imagine 
how the Mixmaster attractor conjecture could possibly be violated. For instance, it seems 
rather absurd that the past attractor consists only of (a subset of) heteroclinic cycles — but 
there exist no proofs. A closely related belief is the following stronger statement. 

Conjecture. For almost all Bianchi type IX solutions 7 the a-limit set otfa) coincides with 
the Mixmaster attractor Aix ■ 

We use the term 'almost all' in a noncommittal way without specifying the measure; recall that 
the word 'generic' already has the well-defined meaning of 'not past asymptotically self-similar'. 
(The usage of the word 'generic' is in accord with [26].) 

Stochastic beliefs 

This subsection is concerned with the (open) question of which role the Mixmaster/Kasner 
map plays in the asymptotic evolution of type IX solutions. The basis for our considerations 
are the results of Section [5] where we discussed the Mixmaster/Kasner map and the stochastic 
aspects of (generic) Kasner sequences. The Mixmaster stochasticity conjecture supposes that 
these stochastic properties carry over to almost every type IX orbit when represented as an 
approximate Kasner sequence. 

Conjecture (Mixmaster stochasticity). The approximate Kasner sequence (tij)igN associated 
with a generic type IX orbit admits a stochastic interpretation in terms of the probability dis- 
tribution associated with the Kasner map, cf. Section^ (This holds with probability one, i.e., 
for almost every generic type IX orbit.) 

The Mixmaster stochasticity conjecture is based on a rather suggestive simple idea: Type IX 
evolution is like trying to follow a path of an (infinite) network of paths while the ground is 
shaking randomly, and where the shaking subsides with time but never stops. A type IX orbit 
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tries to follow a sequence of transitions on the Mixmaster attractor while random errors cause 
the orbit to lose track. Due to the errors, the orbit is incapable of following one particular 
sequence forever; after a finite time, the type IX orbit has deviated too much and it leaves the 
vicinity of the sequence. Although, temporarily, the orbit is contained in a neighborhood of a 
different sequence, it is bound to lose track of that sequence as well eventually. Accordingly 
the type IX orbit is thrown around in the space of Mixmaster sequence with the effect that the 
type IX orbit inherits the stochastic properties of generic Mixmaster sequences. 

In the following we paint a heuristic picture that makes this idea a little more concrete. Con- 
sider a type IX orbit (approximate sequence of transitions) and the associated approximate 
Mixmaster sequence (Pz); £ n. Let (p|°^); e N be the exact Mixmaster sequence with Pq ^ = Po 
and consider the Taub-adapted neighborhood of this sequence associated with some prescribed 
value e > 0. Finite Mixmaster shadowing entails that there exists a finite piece (Pz)z=o,i,...,Li-i 
of the approximate sequence that is contained in the prescribed Taub-adapted neighborhood 
of the exact sequence (p| ^)z 6 n- However, at I = L\, the approximate Mixmaster sequence 
(Pi)igN leaves the prescribed tolerance interval due to the accumulation of errors. Hence, at 
I = Li we reset the system and consider the exact Mixmaster sequence (Pj )j>Li with initial 
data P^ = Pl x - The approximate Mixmaster sequence (P;)z>Li is contained in the Taub- 
adapted neighborhood of the exact sequence (Pp^z^Li up to I = L 2 — I . At I = L 2 another 
readjustment becomes necessary. Iterating this procedure and concatenating the finite pieces 
(Pj we are able to construct a sequence (P;)zeN with the property that the ap- 

proximate sequence (Pz); s n is contained within the Taub-adapted neighborhood (associated 
with e) of (Pi)z 6 N for all I £ N. The sequence (Pz)z s n is a piecewise exact Mixmaster sequence; 
it is exact in intervals [Li, ij+i). The length 5Li = — £; of these intervals grows beyond 
all bounds as i — » oo, because shadowing takes place with an increasing degree of accuracy. 
The error (of the order e) between the approximate sequence (Pz); 6 n and the piecewise ex- 
act sequence (P;)z £ n results from the accumulation of numerous small errors. This obliterates 
the deterministic origin of the problem and generates a 'randomness' that leads to stochastic 
properties. Accordingly we expect that the exact sequences (Pj )j>l< from which (Pz)zgN is 
built constitute a random sample of Mixmaster sequences and thus truly reflect the stochastic 
properties of the Mixmaster map. As a consequence, although the sequence (P;)z e N is only a 
piecewise exact Mixmaster sequence, it possesses the same stochastic properties as a generic 
Mixmaster sequence. Extrapolating this line of reasoning we are able to complete the argu- 
ment and find that the approximate sequence (Pz)z 6 n itself reflects the stochastic properties 
of the Mixmaster map. To emphasize this aspect of stochasticity of (P;)z s n we use the term 
'randomized approximate sequence'. Some comments are in order. 

In our discussion we have assumed implicitly that the approximate sequence (P;)z g N we consider 
can shadow any exact Mixmaster sequence for a finite number of transitions only. This is not 
necessarily the case. A type IX orbit whose a-limit set is one of the heteroclinic cycles (if such 
an orbit exists!) is an obvious counterexample: For every e > there exists L € N such that P; 
is contained in the Taub-adapted e-neighborhood of the Mixmaster sequence associated with 
the cycle. However, we expect that this type of 'infinite shadowing' holds at most for orbits of 
a set of measure zero. 

A more serious limitation of the intuitive picture that we have sketched is illustrated by the 
following related example. Consider a type IX orbit whose a-limit set is the heteroclinic cycle 



depicted in Figure 7(a) (where we note again that the existence of such an orbit is not proven). 
For the associated approximate Kasner sequence (wz)z e N we have ui — ► (f + Vo)/2 as I — ► oo. 
Consider the piecewise exact Kasner sequence (uz); e N that is associated with the piecewise 
exact Mixmaster sequence. Each piece (u[ I '*)/=£;....,£i+i-i is generated from a value = 
[ko; ki, k%, ■ ■ •] = [f ; 1,1,1,1,..., 1, k„\ • • ■]; we have n — > oo as i — > oo. It is evident that 
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these pieces do not form a random sample of Kasner sequences. At best one could conjecture 
(and it is probably safe to do so) that the collection of the remainders [kn , ■ ■ ■] is such a 
random sample. Even so, the stochastic aspect does not carry over to the approximate sequence 
(uz)ieN, since the approximate sequence leaves the neighborhood of each sequence (u| I ' l ); = L i! ... 
before that sequence has entered its stochastic regime (which is characterized by the remainder 
[kn ; ■■•])■ Again we invoke 'genericity' to save the day: We conjecture that, for almost all 
approximate sequences, the sequences (uf^)i>L i are a true random sample of Kasner sequences 
and that the associated stochastic properties indeed carry over to the approximate sequence 
{ui)k=n itself. 

The piecewise exact sequence (P/)/ £ n consists of pieces of length SLi = Li + \ — Li, where Li 
grows beyond all bounds as i — ► oo, because shadowing takes place with an increasing degree 
of accuracy. However, this does not necessarily imply that SLi — * oo as i — > oo. The latter 
property is directly connected with the Kasner map convergence conjecture. 

Conjecture (Kasner map convergence). Almost every generic type IX orbit is associated with 
an approximate Kasner sequence (w;)/gN such that Sui — > as I — > oo, where Sui describes 
the error between iii+i and the value of the Kasner parameter generated from ui by the exact 
Kasner map (|23[) . 

The relevance of the Kasner map for the asymptotic evolution of type IX solutions rests on the 
validity of the Kasner map convergence conjecture. Let us thus give a more detailed discussion 
and present the line of arguments that leads to this conjecture. 

An orbit X(r) (with segments X(n) —* A~(t/+i), I € N) generates a sequence of check points 
• ■ ■ i— ► Pj i— ► P;+i which in turn yields a map ■ • • i— > ui i— > i— >■ • • • . The parameter 

ui+i (associated with P;+i) is generated from ui (associated with Pf) by the Kasner map plus 
an error Sui , see Figure [5] The magnitude of the error Sui depends on the initial data of 
the segment, i.e., on X(ji); equivalently, we may view Sui as a function depending on (i) the 
position of P; on K°, and (ii) the vector X(ti) — P;, which is orthogonal to K° at P;. 

To obtain an estimate for Sui we introduce the order of magnitude of the error, which we denote 
by Aui. We define Aui to be the average of \6ui\ over all vectors A (77) — P; of equal length 
that are orthogonal to K°; alternatively, we use the somewhat 'safer' definition of Am; as the 
maximum of \5ui\. By design, the order of magnitude Aui of the error is a function of two 
variables: (i) the position of P/ on K°, which is invariantly represented by ui, and, instead of 
X{ti) — Vi itself, (ii) the (orthogonal) distance of X(n) from P; (or, equivalently, from K°), i.e., 
||A(t7)-P/|| (= ||A(t ; )-K°||), which we denote by 5 t , see Figure in brief , Aui = (Au)(u ; ,^). 
In the following we investigate the behavior of the function (Aw) (it, 5) under variations of the 
two arguments. 

(i) (Ait) (it, •). Keep u fixed (and assume that u lies in the interval u € (1, 00) so that its image 
under the Kasner map is finite) . Then Am is a function of the distance 5 such that Ait — > as 
<5 — > (which is a simple consequence of the regularity of the dynamical system and continuous 
dependence on initial data). It is reasonable to conjecture that the relation is in fact monotone. 

(ii) (Ait)(-,<5). Keep the distance 5 fixed, so that Ai2 is a function of u. The fundamental 
observation is that Ait becomes unbounded as (a) it — > 00, and (b) u — > 1. Case (a) is due 
to the intricacies of the flow in the vicinity of the non-transversally hyperbolic Taub points 
(where u = 00); recall that it -1 measures the angular distance of the check point from the 
nearest Taub point. For (b) we consider the orbit Q Q — > T Q (for some a) which corresponds 
to u = 1 1—* u = 00. Suppose that P; coincides with Q a (i.e., ui — 1). In general, P;+i will not 
coincide with T a (independently of the choice of 8 — §1); therefore Ai2; = (Au)(ui, 5i) = 00. If 
P; is very close to Q Q , which means that ui is close to 1, a small deviation of A(r) from the 
type II transition emanating from P; will still be small at the end point A(t; +1 ); however, in the 
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vicinity of the Taub point T Q even a small deviation can translate to a large error Sui between 
and (ui — 1) , which in turn results in the asserted blow-up of the order of magnitude of 
the error. The qualitative properties of Au as a function of S and u are depicted in Figure [HI 

Au 




1 *- u 

6 = 1 

Figure 9: We consider a segment of an orbit X(t) with initial check point P associated with the 
Kasner value u. (In accordance with the standard convention for the Mixmaster/Kasner map, 
the terms 'initial' and 'final' refer to a past-directed time.) The quantity 5 denotes the initial 
distance of the orbit from the Kasner circle, which coincides with its distance from P. The 
final check point is not generated from P via the Mixmaster map, but with an error represented 
by 5ii. The figure gives a schematic depiction of the order of magnitude Au of the error in 
dependence on u and 6. Each curve represents (Au)(-,S), i.e., Aii as a function of u for a 
constant value of the initial distance 5; the lower curve is associated with a smaller value of 5, 
the top curve with a larger value. 

Consider a generic type IX orbit X(t) and its representation as an approximate sequence of 
transitions. By Theorem l6.1l each orbit converges to the Mixmaster attractor. As a consequence 
the distance Si = \\X(ti) — K°| (= ||X(t;) — P;||) converges to zero as I — > oo. Therefore, as 
(5; * (Z ^ oo), the sequence of functions (Au)(-,Si) of Figure [5] converges to zero. However, 
this convergence is merely pointwise and not uniform; this behavior is crucial to understand 
the behavior of Sui as I — > oo. 

If a type IX orbit X(t) converges to one of the heteroclinic cycles on the Mixmaster attrac- 
tor, then there exists e such that ttj G (1 + e^e^ 1 ) for all sufficiently large I. On this in- 
terval the function Au of Figure [9] converges to zero uniformly. We therefore obtain that 
Aui = (Au)(ui, Si) — > and thus Sui — > as I — ► oo for this type IX orbit. However, the Mix- 
master attractor conjecture suggests that almost every type IX orbit X(t) has an associated 
approximate Kasner sequence (ui)i£H that is unbounded, hence the general case is not so clear. 
Since (ui)i^ enters the intervals [1, 1 +e) and (e _1 , oo) for any e, Aui = (Au)(u/, Si) (and thus 
Sui) need not necessarily converge to zero. Let us elaborate. 

Let k > be fixed. By the implicit function theorem, the inequality (Au)(u, S) > k is satisfied 
if and only if u < 1 + f K (5) or«> g K ($) for some functions f R , g K , which satisfy f K (S) — > and 
9k(S) — > oo as 5 — > 0. Accordingly, for a given type IX orbit (and its associated approximated 
Kasner sequence u{), we obtain 

Aui = (Au)(ui,Si)>K & ui<l + f K (Si) (35) 

or iii > g K {Si). We call the union of the two intervals (l, l + / re (<5j)J and (g K (8i), oo) the 'hazard 
zone' (associated with epoch number I); obviously, the hazard zone is decreasing with I. Hence, 
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the order of magnitude of the error, Aui, is larger than k at epoch number I if and only if the 
approximate Kasner parameter iii falls into the hazard zone (associated with I). 

In connection with the Kasner map convergence conjecture the question is how often (uj)igN 
enters the hazard zone, i.e., how often (I35|) occurs: finitely many times or infinitely many times? 

Suppose that the Mixmaster stochasticity conjecture is correct. Then almost every sequence 
(u/); 6 n admits a probabilistic interpretation in terms of the probability distribution (|31[) . Ac- 
cordingly, we expect the question to turn into an example of a '0-1 law': There exists two 
alternatives: (i) Type IX orbits (i.e., sequences (ui)i£n) that satisfy (|35p infinitely many times 
are generic (i.e., of measure 1 in the state space); type IX orbits that do not are non-generic 
(i.e., of measure 0); (ii) Type IX orbits that satisfy (|35[) infinitely many times are non-generic; 
orbits that do not are generic. (In both cases, the two sets, the 'generic set' and the 'non- 
generic set', are probably non-empty.) Which of these alternatives is actually realized depends 
on the rate of decay of (5i)i£ft and the decay and growth of f K and g K , respectively, as 6 — ► 0. 
Case (i) is associated with small (subcritical) rates of decay; case (ii) with large (overcritical) 
decay. We conjecture that case (ii) applies, i.e., l|3"5|) is satisfied only finitely many times for 
generic type IX orbits. The reason for this conjecture is the fast rate of convergence of the orbit 
to the Mixmaster attractor which follows from the decay of the functions N1N2N3, cf. (fT0|) . 
and NiN 2 + N 2 N 3 + N 3 N X , cf. [26]; this leads to the expectation that the r.h. side of ([35]) 
represents a rapidly decaying function^ To conclude our line of arguments in favor of the 
Kasner map convergence conjecture, we note that the actual error Sui can be estimated by 
Aui. Therefore, if Am; < k for all I except for a finite set, then also 5ui < k for all I except 
for a finite set. Since k is an arbitrary positive number, the statement of the Kasner map 
convergence conjecture ensues. 

We round off this section with some further remarks on the conjectures. First, we note that one 
might be led to suppose that there could in fact exist type IX solutions for which the statement 
of the Kasner map convergence conjecture is violated, i.e., Sui -/* as I — > 00. In any case, the 
class of these special solutions is at most of measure zero. Second, we note that the statement 
of the Kasner map convergence conjecture, i.e., Siii — > (I — > 00), and the statement 5Li — > 00 
(i — > 00), cf. the discussion on piecewise exact sequences, are expressions of one and the same 
stochastic property. This emphasizes that the two 'stochastic' conjectures tightly intertwine. 
Third, we remark that it is conceivable that 'almost every' in the two conjectures might be 
in fact 'every'. In that case, the non-generic examples of Kasner/Mixmaster sequences, see 
Section would not have any counterparts among the type IX solutions. Of particular interest 
in this context would be to have an answer to question (ii) of the second subsection of Section [7] 
If there do not exist type IX orbits whose a-limit set is one of the heteroclinic cycles, this would 
be a strong indication in favor of 'every' and against 'almost every'. If there exist type IX 
orbits that converge to a cycle, 'almost every' is the best one can aim for in the Mixmaster 
stochasticity conjecture. 

Finally, let us briefly comment on claims of stochastic and chaotic properties of Bianchi type IX 
asymptotic dynamics. Chaotic aspects of the Kasner map and related maps have been studied 
under various aspects [3 [HI \17\ [20]. However, the relevance of these maps for type IX 
asymptotics rests on the two conjectures in this section. If the conjectures turn out to be wrong 
(e.g., if there exists a generic set of solutions that converge to heteroclinic cycles), then none of 
the results on the Kasner map carry over to full type IX dynamics. Numerical investigations, 
see |13) and references therein, reflect Theorem 16.11 and finite shadowing, but it is implausible 
that numerics can possibly shed light on the actual asymptotic limit of type IX solutions. 
Numerical errors are unavoidable and random in nature; these errors generate precisely the type 
of stochasticity the simulation is looking for. Accordingly, numerical studies will necessarily 



12 Additional support for the conjecture comes from toy models that reflect the instability of the Kasner map 
and its consequences. We will come back to this issue in future work. 
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find type IX solutions to exhibit the same stochastic behavior as generic Kasner sequences and 
are thus not suited to explore the validity of the conjectures. 

8 Billiards and billiard beliefs 

In this section we briefly review the Hamiltonian billiard approach to type IX dynamics in the 
spirit of 6, 31] and make contact with the dynamical systems approach. The metric (|2bj) can 
be written as 

4 g = -(det g)N 2 dx° ® dx° + e~ 2bl uj 1 <g> lj 1 + e" 2fc2 Cj 2 <g> Q? + e~ 2b3 u> 3 <g> cl> 3 , (36) 

where a; is an arbitrary time variable that is related to proper time t by a densitized lapse 
function N according to dt — \/det g N dx°; evidently, det g = exp[— 2(6 1 + b 2 + b 3 )]. The 
Hamiltonian for the orthogonal Bianchi type IX perfect fluid models is given by 

U = N ( 3 E G^^P ~ 3R det 9 + 2 P det -9 = . ( 37 ) 

V a,/S / 

where 3 R is three-curvature and p = po(det g)~ <il+w ^ 2 , cf. the remark following Q14p (where we 
recall that dlog(det g) = 6dr). Q a @ is the inverse of the so-called minisuperspace metric Q a p, 

£ Q^vfi := - £ vitf = £ ^ - (J2 va ){j2 w(3 ) ' ( 38 ) 

a, (3 -y^5 a a f3 

i- e -j Sa/3 is a 2 + 1-dimensional Lorentzian metric. The gravitational potential Uq = — 3 R det g 
is given by 

U G = \ (e- 4fel + e- 4 " 2 + e- 4 " 3 - 2 e - 2 ^+ b ^ - 2 e - 2 ^ +b ^ - 2 e - 2 ^ +b1 ^ ; (39) 

where we have set the structure constants n a to one; the potential for the fluid is given by 
U F = 2p det g = 2p exp[-(l - w)(6 1 + b 2 + b 3 )}. For further details see, e.g., [37] or [TO1 
Chapter 10]; note that b a = ~(3 a , which is used in these references; see also [30] and j31j . 

As argued in [HET], b a is expected to be timelike in the asymptotic regime, i.e., J2 a ^b 13 < 
0. Assuming that b a is timelike allows to introduce new metric variables instead of {6 1 , b 2 , b 3 }. 
Defining p 2 = — ^ Q Q a p b a b@ and orthogonal angular metric variables, collectively denoted by 
7, leads to 

53 Gap db a dbP = -dp 2 + p 2 dn 2 h , (40) 

a, 13 

where <if2 2 is the standard metric on hyperbolic space. Making a further change of variables, 

A = logp= ilog(-^g Q/3 &«^) , (41) 

yields 

E ^ ^ = -4 + ^r 2 < = &o _a [--a + A] ■ ( 42 ) 

a, S3 

Choosing the lapse according to N = p 2 leads to a Hamiltonian of the form 

H = 3 [~*X + A\ +P 2 Y. CA exp(-2p«;A(7)) , (43) 

A 
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where u>a(7) denotes certain linear forms of the variables 7 ", i.e., wa{^) — Y]a WApJ^ 3 ', see |31j 
for details. 

The essential point is that one expects that p — * 00 towards the singularity and hence that 
each term p 2 exp[— 2pwA(l)} becomes an infinitely high sharp wall described by an infinite 
step function Qoo{x) that vanishes for x < and is infinite for x > 0. Accordingly, only 
'dominant' terms in the potential are assumed to be of importance for the generic asymptotic 
dynamics, while 'subdominant' terms, i.e., terms whose exponential functions can be obtained 
by multiplying dominant wall terms, are neglected. In the present case there are three dominant 
terms in Uq (which is the minimal set of terms required to define the billiard table) , the three 
exponentials exp(— 46"). Dropping the subdominant terms in the limit p — > 00 leads to an 
asymptotic Hamiltonian of the form 

3 

Woo = \ [-7^ +^]+Yl 0oo(-2 WA ( 7 )) , (44) 

A=l 

where only the three dominant terms appear in the sum. The correspondence between the 
dynamical systems picture and the Hamiltonian picture is easily obtained by noting that the 
Hamiltonian constraint (|37p is proportional to the Gauss constraint @. The dominant terms 
correspond to the terms N 2 , a = 1, 2, 3, in ([9]); the subdominant terms are collected in An- 

The non-trivial dynamics described by (|4~4"|) resides in the variables 7, i.e., in the hyperbolic 
space. It can be described asymptotically as geodesic motion in hyperbolic space constrained by 
the existence of sharp reflective walls, i.e., the asymptotic dynamics is determined by the type IX 
'billiard' given in Figure[Tni Based on the heuristic considerations the limiting Hamiltonian (|44]) 
is believed to describe generic asymptotic dynamics. 



^73 



7i 




Figure 10: The Bianchi type IX billiard: The disc represents hyperbolic space. The asymptotic 
description of a solution is given by geodesic billiard motion inside the triangle, which acts as 
a stationary infinite potential wall, yielding a 'configuration space picture' of the asymptotic 
dynamics. 

The Hamiltonian picture (as represented by the billiard) and the dynamical systems picture (as 
represented by the Mixmaster attractor) are 'dual' to each other. In the Hamiltonian picture 
the asymptotic dynamics is described as the evolution of a point in the billiard. Straight 
lines (geodesies) in hyperbolic space correspond to Kasner states. Wall bounces correspond to 
Bianchi type II solutions; the bounces change the Kasner states according to the Kasner map. 
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Since the billiard picture emphasizes the dynamics of the configuration space variables, one 
may say that the Hamiltonian billiard approach yields a 'configuration space' representations 
of the essential asymptotic dynamics. 

In the dynamical systems state space picture the motion in the Hamiltonian billiard picture 
becomes a 'wall' of fixed points — the Kasner circle K°. The walls in the Hamiltonian billiard 
are translated to motion in the dynamical systems picture — Bianchi type II heteroclinic orbits; 
these type II transitions yield exactly the same rule for changing Kasner states as the wall 
bounces: the Kasner map. Since the variables E Q are proportional to 7r Q , it is natural to 
refer to the projected dynamical systems picture as a 'momentum space' representation of the 
asymptotic dynamics. 

Summing up, 'walls' and 'straight line motion' switch places between the Hamiltonian formu- 
lation and dynamical systems description of asymptotic dynamics, and the two pictures give 
equivalent complementary asymptotic pictures. 

The above heuristic derivation of the limiting Hamiltonian rests on two basic assumptions: It 
is assumed that b a is timelike in the asymptotic regime and that one can drop the subdom- 
inant terms. These assumptions correspond to assuming that £1 and An can be set to zero 
asymptotically, i.e., the procedure precisely assumes Theorem 16. II (The Hamiltonian analysis 
in [351 Chapter 2] uses the same assumptions, and hence the present discussion is of direct 
relevance for that work as well.) That these assumptions are highly non-trivial is indicated by 
the difficulties that the proof of Theorem 16.11 has presented; cf. [26J [27] . An alarming example 
is Bianchi type VIII; in this case the heuristic procedure in this respect leads to exactly the 
same asymptotic results, but so far there exist no proof that f2 and An tend to zero towards 
the singularity for generic solutions. We elaborate on the differences between type VIII and 
type IX in [27]. Moreover, like in the state space picture there exists no proof that all of the 
possible billiard trajectories are of relevance for the asymptotic dynamics of type IX solutions. 
For example, there exists a correspondence between periodic orbits in the dynamical systems 
picture and the Hamiltonian billiard picture, and it is not excluded a priori that solutions are 
forced to one, or several, of these. We emphasize again that all proposed measures of chaos 
that take billiards as the starting point, see [38] and references therein, rely the conjectured 
connection between the Mixmaster map and asymptotic type IX dynamics. 

The above discussion shows that there are non-trivial assumptions and subtle phenomena that 
are being glossed over in heuristic billiard 'derivations.' Nevertheless, we do believe that the 
billiard procedure elegantly uncovers the main generic asymptotic features, and it might be 
fruitful to attempt to combine the dynamical systems and Hamiltonian picture in order to 
prove the Mixmaster conjectures. A tantalizing hint is that in the billiard picture ir\ becomes 
an asymptotic constant of the motion. Rewriting this dimensional constant of the motion in 
terms of the dynamical systems variables yields 



where log V det g — 3(r — To) and To is a constant. 

9 Concluding remarks 

The purpose of this paper is two-fold. On the one hand, we analyze the main known results on 
the asymptotic dynamics of Bianchi type IX vacuum and orthogonal perfect fluid models to- 
wards the initial singularity. The setting for our discussion is the Hubble-normalized dynamical 
systems approach, since this is essentially the set-up which has led to the first rigorous state- 
ments on Bianchi type IX asymptotics . (We choose slightly different variables to emphasize 




(45) 
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the permutation symmetry that underlies the problem.) The main result ('Mixmaster fact') is 
Theorem 16 . 1 1 which is due to Ringstrom [26j : for an alternative proof see [27) . This theorem, in 
conjunction with an analysis of the Mixmaster attractor, leads to a number of further rigorous 
results which we list as consequences [27] . 

On the other hand, we draw a clear line between rigorous results ('facts') and heuristic consid- 
erations ('beliefs'). We make explicit that the implications of Theorem [O] and its consequences 
are rather limited, in particular, the rigorous results do not give any information on the details 
of the oscillatory nature of Mixmaster asymptotics. The mathematical methods required to 
obtain proofs about the actual asymptotic Mixmaster oscillations are yet to be developed; it is 
likely that radically new ideas are needed. This paper provides the infrastructure that might 
yield the basis for further developments. Our framework enables us to transform vague beliefs 
to a number of specific conjectures that describe the expected 'complete picture' of Bianchi 
type IX asymptotics. The arguments we give in their favor are based on an in-depth analysis 
of the Mixmaster map and its stochastic aspects in combination with the dynamical systems 
concept of shadowing. 

We conclude with a few pertinent comments. First, as elaborated in [27) there exists no cor- 
responding theorem to Theorem 16.11 in other oscillatory Bianchi models such as Bianchi type 
VI_i/g or VIII; this suggests that the situation in the general inhomogeneous case is even more 
complicated than expected. Furthermore, numerical studies are incapable of shedding light on 
the asymptotic limit. This is mainly due to the accumulation of inevitable random numerical 
errors that make it a priori impossible to track a particular type IX orbit. Finally, although the 
Hamiltonian methods are a formidable heuristic tool, so far this approach has not yielded any 
proofs about asymptotics. Nevertheless, it might prove to be beneficial to explore the possible 
synergies between dynamical systems and Hamiltonian methods. 

In this paper and in [27[ 130) we have encountered remarkable subtleties as regards the asymp- 
totic dynamics of oscillatory singularities; this emphasizes the importance of a clear distinction 
between facts and beliefs. 
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